Exact spontaneous plaquette ground states for high-spin ladder models 
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We study the exchange physics in high spin Mott insulating systems with S — 3/2 which is real- 
izable in ultracold atomic systems. The high symmetry of 5*0(5) or 57/(4) therein renders stronger 
quantum fluctuations than the usual spin- 1/2 systems. A spontaneous plaquette ground state with- 
out any site and bond spin orders is rigorously proved in a ladder spin-3/2 model, whose topological 
excitations exhibit fractionalization behavior. The generalization to the SU(N) plaquette state is 
also investigated. 
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There has been considerable interest recently in high 
spin systems with enlarged symmetry among both 
condensed-matter and atomic physic s 1 ^ 3 : 4 : 5 : 6 . With the 
rapid progress in the field of ultra-cold atomic physics^, 
optical traps and lattices open up new possibilities of 
simulating and manipulating high spin physics experi- 
mentally in a controlled way. For example, the high spin 
fermionic 40 K gas has been produced in one-dimensional 
optical lattices 8 . On the other hand, as a paradigm in the 
low-dimensional magnetic systems, the spin chains and 
ladders have stimulated intense investigation recently of 
how to simulate them using cold atoms. Several schemes 
have been proposed to implement these quantum spin 
models in optical lattices^uLiS. So far most of these pro- 
posals are still concentrated in the spin- 1/2 and spin-1 
systems. Although most high spin systems only exhibit 
the usual spin SU(2) symmetry, a generic high symme- 
try of 50(5), or isomorphic Sp(4), has been rigorously 
proven m spin 3/2 systems 1 ^. They may be realized 
with several candidate atoms, such as 132 Cs, 9 Be, 135 Ba, 
137 Ba. This symmetry sets up a framework to unify 
many seemingly unrelated properties of Fermi liquid the- 
ory, Cooper pairs, and magnetic structures in such sys- 
tems. Conventionally, a high spin system is assumed in 
the large- 5 limit, and thus is considered as more classical- 
like than its low spin counterpart. However, due to their 
high symmetry, the spin-3/2 systems are actually in the 
large- N limit. As a result, the quantum fluctuations are 
stronger than the usual spin 1/2 systems, which results 
in many novel properties. 

In this work, a spin-3/2 exchange model with intrinsic 
50(5) symmetry is proposed by us for the first time, and 
it includes the SU (4) model as its special case. We then 
study a class of solvable spin-3/2 ladder models which 
exhibit the exact spontaneous plaquette ground states 
without any site and bond spin orders, and the topolog- 
ical excitations are studied. To the best of our knowl- 
edge, the existence of a plaquette phase has never been 
exactly proved before, therefore our results also provide a 
firm ground for understanding the plaquette phasa 4 * 1 ^* 1 ^. 
The generalization of our theory to arbitrary spin systems 
with SU (N) symmetry is also discussed. 

We first derive the general 50(5) exchange model from 



a spin-3/2 Hubbard model in the strong repulsive inter- 
action limit Uq,U 2 3> t, where C/0,2 are the on-site Hub- 
bard repulsion in the singlet and quintet channels, respec- 
tively. The projection perturbation theory is employed 
to study the low energy exchange process through vir- 
tual hopping at quarter-filling, i.e. one particle per site. 
For two neighboring sites, the total spin can be Stot = 
0,1,2,3. Exchange energies in the singlet and quin- 
tet channels are Jo = 4t 2 /C/ , J 2 = At 2 /U 2 respectively. 
These two channels also form 5p(4)'s singlet and quin- 
tet. On the other hand, no exchange energies exist in the 
triplet and septet channels, which together form a 10D 
representation of the Sp(A) group. The effective Hamil- 
tonian can be expressed through the bond projection op- 
erators in the singlet channel Qo(*j) and quintet chan- 
nel Q 2 (ij) as H ex = Y^{i,j) h H = {JnQoihj) + 
J 2 Q 2 (i, or in an explicitly Sp(A) invariant form as 

= E ij IciA^(i)A^(j)+c 2 A^(i)A^(j)] ~ c 3 } 
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4 , C 2 = U "Y J ° aild C 3 

A 11 = 2L ab = cllfpcp (71 = 1, • • • , 10, 1 < a < b < 5) 
and A'' 2 = 2n a = 4 T a a0 cp ( 72 = 11, • • • , 15, 1 < a < 5) 
with the Dirac T matrices which could be found in Refi 
The n a operators transform as a 5- vector under the 5p(4) 
group and the L a \, operators form the 10 generators of 
the 50(5) group. 

It is obvious that an 5t/(4) symmetry appears at Jo = 
J 2- Then Eq. (JIJ reduces into the 5£/(4) symmetric form 
with the fundamental representations on every site where 
L a h and n a together (or A 1 with 7 = 1, • • ■ , 15) form the 
15 generators of the SU (4) group. In this case, the spin 
3/2 Hamiltonian is equivalent to the well-known SU(4) 
spin-orbital model up to a constant ferm^i 1 .. Equation 
JU can be expressed by the usual SU(2) spin operators 
with bi-quadratic and bi-cubic terms as 
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a (5, • Sj) + b(Si- Sj) 2 + c {Si ■ 5j) 3 , (2) 
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FIG. 1: (a)The spin ladder model. Plaquettes in bold lines 
represent the plaquette singlets. (b)The basic block to con- 
struct the spin ladder model. 

C = ll( J + J 2)- 

The magnetic structure of the one-dimensional spin- 
3/2 model |J2J can be analyzed by bosonization meth- 
ods combined with renormalization group methods via 
studying the corresponding Hubbard model at quarter- 
filling 13 . A spin gap opens with the appearance of the 
spin Peierls order at XJi > Uq > 0, while it varnishes at 
Uq > U2 > O^, where /Jo, 2 are the on-site repulsions in 
the spin singlet and quintet channels respectively. 

Now we construct a spin-3/2 ladder model exhibiting 
an exact plaquette ground state, which is an SU(4) gen- 
eralization of the valence bond state (VBS) as shown in 
Fig 1. Different from the usual SU{2) case where ev- 
ery bond is an SU{2) singlet, here at least four sites, i.e. 
a plaquette, are needed to form an 5(7(4) singlet, thus 
there is only plaquette order without any bond order in 
this state. As shown in Fig. 1(a), the model is described 
by the Hamiltonian 

Hlad = J± / ] fyi,l),(i,2) + J /J fyi,q),(i+l,a) 
i i,a 
2 

+ ^2^2 J i [h(i,l),(i+6,2) + fyi,2),(i+(5,l)] 
i 5=1 

(3) 

i,a 

where the site is labelled by its rung num- 
ber i and chain index a — 1,2 and h{j = 
\ [ciA"* 1 (i)A 71 (j) + c 2 A~> 2 {i)A 12 (j)] . Here we take 
J± , J, Jy , J y and J' as positive constants in unit of 
Jo and set J = 1. The regular railway ladder model 
corresponds to J 1 = Jy = J 1 = 0. To gain some 
intuition about the plaquette state, we first consider a 
4-site system with diagonal exchanges (a tetrahedron). 
For such a simple system with SU (4) symmetry, the GS 
is a plaquette singled defined as 

SU4 (1234)= -J^f^ |1 P 2,3 7 4 5 ) 

where T^^s is an antisymmetric tensor and 
/x, v,7, S — ±|,±2< Such an SU(4) singlet is rota- 
tionally invariant under any of the fifteen generators A 1 



of an SU(4) group with i = 1, ■ • ■ ,15. For a tetrahedron 
with SO(5) symmetry, the SU(4) singlet is no doubt 
an eigenstate with the eigen-energy = — 5ci — 5C2. 
Furthermore, it is easily verified that the 5J7(4) singlet 
is the ground state of the four-site SO (5) exchange 
model for C2 > 0. 

We now focus on the SU(4) ladder described by Eq. 
© with c 2 = ci. It turns out that the exact ground 
state of the SU(4) ladder model is a doubly degenerate 
singlet provided the relation 

J± = l-r= |/i = 3J*=3J' (4) 

is fulfilled. The degenerate GSs are composed of products 
of nearest-neighboring plaquette singlets. Explicitly, for 
a ladder with length M, the two degenerate ground states 
are given by 

M/2 

\Pi) = n su 4[(2i,l),(2j,2),(2 l + l,l),(2 l + l,2)] 

i=l 
M/2 

|P 2 ) = n su 4[(2z-l,l),(2*-l,2),(2z,l),(2z,2)] 

i=l 

with the corresponding ground energy 

E = -M^c x , (5) 

where periodic boundary condition is assumed M+l = 1. 

The proof of the above conclusion can be understood 
through two steps. First, one observes that the plaquette 
singlet product is no doubt an eigenstate of the global 
Hamiltonian because any generator of the ST/ (4) aside 
from the plaquette acting on the plaquette singlet is zero, 
i.e. 

A 5 [ A J + A l + A a + A J\ 5u 4 (1234) |5i/) = 

where 7 = 1, • • • , 15. Secondly, we prove that such an 
eigenstate is the ground state of the global Hamiltonian. 
In order to prove it, we utilize the Rayleigh-Ritz varia- 
tional principle 

M M 

(*Ch(50l*) > E 9-s. >£>»•'• (6) 

i=l i=l 

which implies that an eigenstate |\E') is the ground state 
of a global Hamiltonian if it is simultaneously the ground 
state of each local sub-Hamiltonianii Here E g s . is the 
ground state energy of the global Hamiltonian which 
is represented as a sum of M sub-Hamiltonians, say, 
H = ^2iLih(Bi), and e g . s .(Bi) represents the ground 
state energy of a sub-Hamiltonian h(Bi). 

Now we apply the above general principle to study our 
model given by Eq. J3J . Explicitly, as long as Eq. (0} holds 
true, we can decompose Eq.Q as 

M 

H = Y J h{B l ), (7) 

t=i 
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where h(Bi) = J' h-ij denotes the Hamiltonian of a 
six-site block as shown in Fig. 1(b) and (ij) represents all 
the available bonds in the block of Pj. For convenience, 
we use the rung index on the left side of a six-site block to 
label the block. For such a six-site cluster with 15 equiv- 
alent bonds, the local Hamiltonian can be represented as 
a Casimir operator and the representation with smallest 
Casimir corresponds to the Young diagram [2 2 1 2 ]. It 
follows that eigenstates given by 

\Bi) x = «u4[(i,l),(i s 2),(i + l,l),(i + l,2)]®dn-2 
\Bi) 2 = d i ®«t 4 [(t + l,l),(* + l,2),(i + 2,l),(i + 2,2); 

are two of the degenerate ground states of the six-site 
block Bi with the ground energy e g {Bi) — — J'(^ci + 
|ci), where di represents a dimer on the ith rung 

di = [(*, 1), (i, 2)] = -^V | (i, 1)^,2),,), 

with T^v denoting an antisymmetric tensor. A global 
eigenstate can be constructed by a combination of local 
eigenstates \B i ) 1 or \Bi) 2 because the dimers on the side 
of the block are free in the sense that they can form pla- 
quette singlets with other sites belonging to the neighbor- 
ing blocks. Now it is obvious that the eigenstates \P\) 
and \P 2 ) are the ground state of each sub-Hamiltonian 
h(Bi), and therefore the degenerate GSs of the global 
Hamiltonian with the GS energy given by Eq.JSJ). In fact, 
the constraint relation Eq. (0} can be further released to 

J± ~ \ J = l J * = 3J ' = 3J '' (8) 

which means that the degenerate plaquette products are 
the GSs even in the strong rung limit J± 3> J. Further- 
more, we note that \P{) and \P 2 ) are the eigenstate of 
the SU(4) ladder model even for arbitrary J±, but not 
necessarily the GS. 

This plaquette-product state is a spin gapped state 
with only short-range spin correlations. The two GSs are 
spontaneously tetramerized and thus break translational 
symmetry. The elementary excitation above the GS is 
produced by breaking a plaquette singlet with a finite 
energy cost, thus leading to an energy gap. Two kinds of 
excitations are possible, either a magnon-like excitation 
or a spinon-like excitation. The gap size of a magnon- 
like excitation is approximately equivalent to the energy 
spacing level between the GS and the first excited state of 
a local tetrahedron. For the model JSJ, the first two ex- 
cited states above the plaquette singlet are represented by 
the Young diagrams [2 : 1 2 ] and [2 2 ] and the correspond- 
ing gap sizes are A m = 4ci and A s = 6ci respectively. 
The spinon-like excitation is made up of a pair of rung 
dimers which propagate along the leg to further lower 
the energy, and thus the excitation spectrum is a two- 
particle continuum. The propagating dimer pairs behave 
like domain- wall solitons (kink and anti-kink) connecting 
two spontaneously tetramerized ground state. We repre- 
sent an excited state with a kink at site 2m — 1 and an 




FIG. 2: Two-fold degenerate ground state of the SU{A) 
diamond-chain model. The Heisenberg exchange is denned 
on every linked bond with the equal magnitude. 

anti-kink at site 2n as $ (m, n) , so the corresponding 
momentum-space wavefunction is 

V(ki,k 2 ) = e l(2m - 1)kl+l2nk ^ (m,n) . 

l<m<n<M 

The excited energy can be calculated directly by using 
the above variational wavefuctionM. For a spontaneously 
tetramerized ladder system, it is reasonable to assume 
that the kink and antikink are well separated because 
there exists no intrinsic mechanism of binding them to- 
gether to form a bound state, therefore we treat kink and 
antikink separately. Since the state \P (m) is not orthog- 
onal with the inner product given by (\& (m')\ |\P (rn)} = 




, thus \& (fci , k 2 ) has a nontrivial norm. After 



considerable but straightforward algebra, we get 
37 6 

e{k 1 ,k 2 ) = — A - — (cos 2k x + cos 2k 2 ) A, 
35 35 

where A = 5ci is the energy level spacing between the GS 
and the excited state where a local plaquette is broken 
into two rung singlets. Since the SU(4) spin-3/2 model 
can be mapped into the spin-orbital model, our above 
results can be directly applied to the corresponding spin- 
orbital ladder model. For an SU(4) spin-orbital railway- 
ladder model, we noticed that it has similar degenerate 
GSs composed of SU(4) singlet from the work of Bossche 
et al. based on the exact diagonalization method and 
analytical analysis in strong coupling limits^. 

For the SO (5) ladder © with c 2 > 0, one might expect 
that the degenerate plaquette states are the GSs in the 
parameter regime (jHJ) like in the case of the SU(4) ladder. 
Unfortunately, one can check that those plaquette states 
are no longer the eigenstates and therefore not the exact 
GSs. Nevertheless, we can prove that |Pi) and \P 2 ) are 
degenerate eigenstates of the global Hamiltonian if J± = 

J J y-^ 2 J 2 t/ . 

Next we turn to a solvable diamond-chain model as 
shown in Fig. 2, where the spin-3/2 Heisenberg model is 
defined on the linked bonds with the same exchange en- 
ergy H = J~Y^Uj) hij- For an SU(4) model, the Hamilto- 
nian can be written as a sum of the Casimir of the to- 
tal spin in each five-site cluster. Thus we conclude that 
the GSs are doubly degenerate plaquette singlets because 
among the possible representations composed by 5-sites, 



4 




FIG. 3: Excitations in the 5(7(3) model as the nearest- 
neighbor Heisenberg model defined on the ID corner-shared 
tetrahedrons. 

the above ground state configurations provide the small- 
est Casimir. By the same reasoning as above, this con- 
clusion holds true for the SO(5) case with C2 > 0. The 
excitation is composed of a spinon and a three-site bound 
state as a result of a plaquette singlet breaking up to a 
1+3 pattern, which bears some similarity to the above 
spin ladder model where symmetric spinon-like excita- 
tions correspond to a plaquette singlet breaking up to a 
2+2 pattern. 

The exact plaquette GS discussed above can be 
straightforwardly generalized to the SU(N) case. The 
SU{2) case corresponds to the celebrated spin- 1/2 
Majumdar-Ghosh modeli^*^ The ST/ (3) case as shown 
in Fig. |3| is particularly interesting. It is defined in a 
ID correspondence of the pyrochlore lattice where only 
the nearest neighbour bond interaction is involved. The 
Hamiltonian is 

ff = J^^^(^(j), (9) 

(«> 7=1 

where A's are the eight SC/(3) generators in the funda- 
mental representations, and (ij) means the sum over the 
nearest neighbors. It can also be written as the spin-1 
representation as H = JY^(a) ^(^)'S(j)+a(S(i)-S(j)) 2 , 



where the SU(3) point is located at a = 1. The Hamilto- 
nian can be written as a sum of the Casimir of the total 
spin in each tetrahedra. We know that the representa- 
tions with the smallest Casimir made out of four sites 
in the fundamental representations is three dimensional. 
The excitations also have a gap and are fractionalized 
as in the ID polymer systems. The three-site singlet 
(quark) is broken up to a 1+2 pattern as one monomer in 
the fundamental representation, and the other two form 
an anti-fundamental (anti-quark) representation. Thus 
the quark and anti-quark pair states are 3*3 fold degen- 
erate. The singlet site monomer can hop around in the 
two faces without breaking more singlets. The detail will 
be the subject of a future publication. 

In summary, we derived and studied the effective 
SO(5) spin-3/2 exchange model as well as the spin ladder 
models with exact two-fold degenerate plaquette GS. Our 
results indicate the formation of spontaneously tetramer- 
ized GS for a translational invariant spin ladder system. 
Quantitative results for the elementary excitation spec- 
trum of the ST/(4) spin-orbital ladder are also obtained. 
Due to the existence of an intrinsic SO(5) symmetry, we 
expect that the plaquette phase of the spin-3/2 models 
can be observed in optical lattices in the future exper- 
iment. The generalization of our theory to the system 
with SU (3) or even SU (N) symmetry is also addressed. 
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